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Abstract: Material and its processing effect the relative structure contains small flaws whoa* 
distribution. These may vary from nonmetallic inclusions and micro voids to weld defects, gji 
quench cracks, surface laps, etc. The objective of a Fracture Mechanics analysis is to det 
small flaws will grow into large enough cracks to cause the component to fail cata^coprHfally. By using 
fracture mechanics theories i.e. LEFM approach, applied to two dimensional objea^UjA as plates and 
shells etc. The aim of this paper is to estimate numerically the beneficial effects? \ajHate with hole and 
bounded single and double side patches are symmetrically bonded on the ton^fi^bottom on surface for 
repairing crack using the finite element method. An observed the change in stfaeaintensity factor along the 
crack length by using FRANC-2D/L simulation software. The obtained iea^ts^now the reduction of stress 
intensity factors. The effects of the adhesive, patch and plate prop^ft^l^yi the rate of reduction of the 
stress intensity factors are highlighted. 
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1. Introi 

The determination of fracture mechanics parameteA^Jich as Stress Intensity Factors (SIF) (Kj, Kn, Km) 
which corresponds to three basic modes of fractMS^-Integral a path-independent line that measure the 
strength of the singular stresses and strains neara**rack tip. Energy release rate (G) which represents the 
amount of work associated with a Crack Q]£rHpg Displacement (COD) or Closure. By comparing these 
parameters with critical values one can eiuH^re the brittle failure state of structures. The crack analysis is 
found on three different following met™dJ[2]. Which are displacement correlation methods, virtual crack 
extension methods, and modified cl^^^closer method. The FRANC-2D/L simulation software is based on 
Finite element method (FEM) ara^bouhdary element method (BEM) are the most widely use techniques for 
evaluating stress intensity facra^SlF). The most important region in modeling the fracture is around the 
crack. While the domain i^Cgied we are using crack tip elements with singularity [16, 17], those elements 
remove the nodal singitfarky at the crack tip. Displacement correlation method was employed. To 
determine stress inter/frwSctors (SIF) 

The main objecntl^f/ihis paper is the stress intensity factors are calculated from displacement correlation 
method. 

1. Displacement Correlation Method 

Bfoplawement correlation is a direct approach which the simplest and historically one of the first technique 
usel^ro estimate SIFs from FEM results [3]. In its simplest form, the finite element displacements for one 
point in the same mesh are substituted directly into the analytical expressions for near tip displacements, 
after subtracting the displacements of the crack tip [1] . 

The point is to be selected as a node on the crack face where the displacements will be greatest, and thus 
the relative error in the displacements is expected to be smallest. The configuration for this simple 
approach is shown in fig 1. 
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Fig.i. Possible sample point locations for displacement correlation. 



The expressions for the SIF's using plane strain assumptions are: 
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The same expressions can be used for plane stress assumptions o^^Nsweplaced with: 

3 = 9/ 



7(l + 3) 

J_ 

Crack tips produce a ^ singularity. The strafe fields near a crack tip of an isotropic linear elastic material 

7" ^ 



can be expressed as a product of 



f^Ofu i 



i of 9 with a scaling factor K: 




Fig .2. Stress field near crack tip of an isotropic linear elastic material. 
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2. Problem Outline 

The objective of this work is to illustrate the process of finite element program to compute the stress 
intensity factor histories accurately so that they can be used for fatigue crack growth (FCGR) antWfe 
predictions. 



Mode I (Tension. Opening) 




Fig.4. Opening mode of deformation. 

Fig.5 shows a simple plate with hole under symmetric 1 
(considering the right half). The traction load is applied* 
initiated at the location of high effective stresses^ 
FRANC 2D/L (Cornell Fracture Group, 2002) [6] i 




'Tttttl 

Fig.5. |j£Je width plate with hole 



Let us consider half portion of the Figure 
ktne top and bottom edges of the plate. A crack is 
lows it to propagate in a straight line manner, 
the finite element solver. 



:(t)«cg 



4? 



The input parameters are: 
Young's modulus (E) = 10.06E6 N/rn 
Poisson's ratio ( i9 ) = 0.25 
Thicloiess of the plate (t) ^ 
Radius of the hole (r) ^ 
Width of the pla^^^^4mm 
Height of thitfQrV = 8mm 
Assumat&h#are: Material is linearly isotropic, 
T^^^pblem is modeled as a Plane stress condition. 

3. Approach 

As discussed, the problem is reduced based on symmetry, the plate is modeled by using CASCA [7], and 
then a converged mesh is made with 8-noded quadrilateral element as shown in fig7. A traction load of 
10000N is applied along the top and bottom edges of the plate. This is the condition for the simple plate 
with hole without addition of any patches and the results are as shown in tablei, which will be discussed in 
results and discussions. 
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Fig.6.Mesh for plate 



Fig.7. Patch mesh 



The mesh for patch after adding patch is as shown in fig 8. The different layers patch 
adhesive are as shown in fig 9 and 10. The present case is the case for the plate with 
In this case on patch is added at the bottom of plate. 



For double patch, the patches are added at the top and bottom sides of the plai^k- 



plate after adding 
*th single patch. 




Fig.8. Plate Layer 




Fig.9. Patch 



Fig.11. Adhesive elements 
s and Discussions: 

The following table shows the vaduesa£5lF obtained for plate with hole without any patches 
Tablejo^ws the results for plate with hole with single patch. 



<5 



si no 


crack length (mm) 


N/mm 2 -nrai' 2 




0.5 


0 


2 


0.6 


14660 


3 


0.8 


19300 


4 




22710 


5 


1.2 


27170 


6 


i-4 


33260 



Table2. Shown below shows the results for plate with hole with single patch 



si no 


crack length (mm) 


N/mm 2 -nrai' 2 




0.5 


4246 
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0.6 


4155 


3 


0.8 


4247 


4 




4256 


5 


1.2 


4372 



Table3. SIF values for plate with hole with double patch 




The following graph shows the variation of SIF along the crack lengj 



ie above three cases discussed. 



Fig6: Varjtt&m^& tress intensity factor V/s crack length 
The following graph is obtained Aom experimental observations. 

J> IS 



<5 




Crack Length (inches) 

Fig 13: Variations of the stress intensity factors (experimental results) 
Conclusions 



In this work The problem is considered under mode-I conditions and K] the stress intensity factors are 
calculated from displacement correlation method .An observed the change in stress intensity factor along 
the crack length by using FRANC-2D/L simulation software for a plate with hole and bounded single and 
double side patches are symmetrically bonded on the top and bottom on surface for repairing crack using 
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the finite element method. The increasing the number of patch layers reduces the stress intensity factors 
which mean a safe design. The results are plotted by varying crack length are compared with the 
experimental graphs .The obtained results are found to be satisfactory 
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